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ABSTRACT 

This review identifies a number of the principal research advancements that have 
occurred over the last five years in the study of electromagnetic (EM) waves in the 
Earth’s inner magnetosphere. The observations used in this study are from the plasma 
wave instruments and radio sounders on Cluster, IMAGE, Getotail, Wind, Polar, 
Interball, and others. The data from passive plasma wave instruments have led to a 
number of advances such as: determining the origin and importance of whistler mode 
waves in the plasmasphere, discovery of the source of kilometric continuum radiation, 
mapping AKR source regions with “pinpoint” accuracy, and correlating the AKR source 
location with dipole tilt angle. Active magnetospheric wave experiments have shown 
that long range ducted and direct echoes can be used to obtain the density distribution of 
electrons in the polar cap and along plasmaspheric field lines, providing key information 
on plasmaspheric filling rates and polar cap outflows. 


1. INTRODUCTION 

Wi thin the magnetosphere there are more then 40 plasma wave emissions that have been 
classified (see the excellent review by Shawhan,\919). It is well known that space 
plasmas can support a variety of EM wave modes. To illustrate. Figure 1 shows various 
propagation modes for a point in the inner magnetosphere (at radial distances less than 
geosynchronous) in which the local electron gyrofrequency is less than the local electron 
plasma frequency (f g < f p ; i.e. the plasmasphere). As the frequency increases (from Hz to 
kHz) plasma waves can propagate in the whistler, Z, and the L-O and R-X modes. The 
frequency range of each mode has been shaded to easily distinguish between them. 

These modes have lower and, sometimes, upper frequency limits to their propagation that 
are called cutoffs or resonances as described by Stix [1992]. The left hand ordinary (L-O) 
and the right-hand extraordinary (R-X) modes are called the free space modes (using the 
terminology of Stix, 1992) since they only have lower frequency cutoffs (fR for R-X 
mode and f p for L-O mode as shown in Figure 1) and at sufficiently high frequencies can 
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propagate into the solar wind once they are generated in the magnetosphere. The Z and 
whistler mode emissions are referred to as trapped emissions since they have well defined 
upper frequency limits that confine these emissions to specific regions in space. 

Due to space limitations it is impossible to review all the progress made in understanding 
magnetospheric processes from EM waves, so the authors will concentrate on a few 
selected topics that have undergone rapid advances in the last five years in both trapped 
and free space portions of the spectrum. Recent missions such as Cluster and IMAGE, 
along with older missions such as Geotail, Wind, Polar, and Interball that routinely orbit 
the Earth, have provided a wealth of new data for analysis. In addition to passive 
measurements, the long range EM sounder capability on IMAGE and the short-range 
sounders on the four Cluster spacecraft are producing unique types of data. The purpose 
of this paper is to review a number of important advances in our understanding of inner 
magnetospheric processes from passive and active wave measurements. The paper will 
first concentrate on passive measurements consisting of trapped emissions (Section 2) 
and escaping emissions (Section 3) and then on active magnetospheric sounder 
measurements (Section 4). 

2. TRAPPED EMISSIONS IN THE PLASMASPHERE 

The index of refraction in the whistler mode is such that some of these EM waves have a 
natural tendency to travel along a magnetic field. In the plasmasphere, whistler mode 
waves make repeated journeys from northern to southern hemisphere and vice versa. 

Since the upper frequency of the whistler mode is the local f p or f g (see Figure 1) 
whichever is less, these waves are largely confined or trapped in the higher density 
plasmasphere region around the Earth. 

There are four basic whistler mode waves generated in the plasmasphere as summarized 
in Figure 2 as a function of their frequency range. From 10-50 kHz, Earth-based 
transmitters dominate the spectrum. The authors have identified over 100 stations 
generating emissions at more than 60 discrete frequencies, and there are quite likely 
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many more. Chorus emissions, observed just outside the plasmasphere typically in the 
frequency range from 0.3 - 12 kHz, can be found extending from the magnetic equator to 
mid-latitudes, and generally intensify during storm periods. Plasmaspheric hiss is 
observed in the plasmasphere at frequencies typically between 0.3 - 3 kHz and is 
believed to be responsible for electron precipitation in the slot region (L-2.5) between the 
inner and outer radiation belts. Electromagnetic (EM) equatorial waves are typically 
found from a few Hz to about 300 Hz and are confined to the magnetic equator (typically 
within 1 R e ) within the plasmasphere. Lightning whistlers are also observed in the 
plasmasphere (not shown in Figure 2) and they have the frequency range of from about 0 
to 10 kHz. 

The average magnetic wave field spectral density measurements, as a function of 
frequency and L value, from the Plasma Wave Instrument (PWI) on Dynamics Explorer - 
1 (hereafter DE) is shown in Figure 3 [adapted from Andrt et al , 2002]. The top and 
bottom panels of Figure 3 are from quiet and disturbed conditions respectively, use the 
same intensity levels, and show the dramatic increase in intensity of the whistler mode 
wave spectrum in the plasmasphere (except for the transmitters) during disturbed 
conditions. The longitudinal distributions, where the whistler mode waves are the most 
intense in the plasmasphere, are shown in Figure 4. Most of the whistler mode wave 
research in the last 5 years has dealt with various aspects of plasmaspheric hiss and 
chorus. Resent results of plasmaspheric hiss research will be discussed in this review 
since chorus is well covered by another paper in this monograph [see Meredith et al, 
2004]. 

2.1 PLASMASPHERIC HISS 

Although found essentially everywhere in the plasmasphere at some intensity level, 
plasmaspheric hiss is most intense throughout the local afternoon sector (see Figure 4) 
and on L shells which contain the slot region in the electron radiation belts. The origin of 
this emission is still somewhat controversial as either generated by gyroresonance or by 
lightning. Research in this area over the last 5 years supports the idea that plasmaspheric 
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hiss is an accumulation of many whistlers trapped in the plasmasphere, but the relative 
contribution from different sources is still unknown. 

The classic theoretical work by Kennel and Petschek [1966] held that plasmaspheric hiss 
receives its energy from a gyroresonance interaction with inner radiation belt electrons 
near the magnetic equator causing the electrons to change pitch angle and precipitate, 
generating the electron slot region between the inner and outer belts. The theory states 
that the amplification of the emission is maximized when the wave normal angle (WNA), 
which is the angle between the wave’s phase velocity (k vector) and the local magnetic 
field, is near 0°. The amplification expected, however, was only a few dB and well below 
the observed intensities of hiss (see Figure 3). Thorne et al [1979] suggested that 
plasmaspheric hiss would only grow in intensity from the background thermal noise to its 
observed intensity from gyroresonance acceleration as the whistler mode wave returned 
through the equator repeatedly. This idea would also require the waves to maintain a 
near parallel WNA. More recently, Able and Thorne [1998] have suggested that this 
would happen if the hiss would suffer reflection from the plasmapause at mid or higher 
latitudes during its inter-hemisphere trajectory. But, this idea is not universally held 
since Storey et al [1991] observed hiss propagating at large WNA even in the equatorial 
region. 

Using the wave distribution function technique developed by Storey and Lefeuvre [1979], 
Santolik et al. [2001] found that most of the observed waves on L shells above 3 had 
WNA nearly parallel to the Earth’s magnetic field, while others had oblique WNAs. In 
addition, these authors found evidence of wave amplification near the geomagnetic 
equator supporting the theories by Kennel and Petschek [1966] and Thorne et al., [1979]. 

Bortnik et al [2003] have found that the lightning generated whistlers tend to settle on 
preferred L-shells in the plasmasphere with the lower frequency components settling at 
higher L-shells and higher frequency components oh lower L-shells. These whistlers are 
referred to as magnetospherically reflected since they do not bounce off the plasmapause 
but are internally reflected within the plasmasphere. By combining the lightning. 
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whistler lifetimes with the power spectral density of lightning Bortnik et al. [2003] 
showed a clear maximum in wave energy in the slot region at about an L of 2. 

Rodger et al [2003] extended these results and determined, from modeling calculations, 
that electrons in the ~50 to 150 keV range precipitate out of the slot region (L = 2 -2.4) 
through gryoresonance interaction with lightning generated whistlers. For electron 
energies above this range, Rodger et al [2003] believes that ground-based transmitters 
and plasmaspheric hiss should dominate over all other loss processes. 

3. CONTINUUM RADIATION 

The Earth’s nonthermal continuum radiation is observed over a veiy broad frequency 
range from as low as 5 kHz [Gurnett, 1975] to as high as 800 kHz [Hashimoto et al 
1999] and is generated in the free space L-O mode above the local f p from sources at or 
near the plasmapause [Gurnett et al., 1988]. In the last 5 years, research in continuum 
radiation has focused on improving our understanding of the source location, emission 
cone characteristics, and detailed spectral measurements. Much of what has emerged 
from these studies in terms of source location is summarized in Figure 5. The lower 
frequency trapped and escaping continuum is generated in the pre-noon section (also 
called normal continuum), the continuum enhancement is generated in the morning 
sector, and the kilometric continuum is generated in deep plasmaspheric notch structures. 
These emissions will be discussed below. 

Continuum radiation has been shown to be closely associated with strong narrow-band 
electrostatic emissions at the plasmapause at the magnetic equator [Kurth et al. , 1979]. 
These strong electrostatic bands occurs at frequencies where the frequency of the 
electrostatic upper hybrid resonance (fuhr) is equal to the frequency of the electrostatic 
(n+l/2)f g resonance. It is believed that continuum radiation is generated by a mode 
conversion process, converting electrostatic Z mode emissions to the electromagnetic L- 
O mode, over the entire frequency range of the emission (see Figure 1). The two classes 
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of mode conversion mechanisms that have been considered are linear [e.g., Jones , 1976 
and Sudden , 1980] and non-linear [e.g., Melrose , 1981; Fung and Papadopoulos , 1987]. 

3.1 NORMAL CONTINUUM RADIATION 

At frequencies less than the magnetopause plasma frequency (~ 30 kHz), the continuum 
radiation has been referred to as the "trapped" component [Gurnett and Shaw , 1973] 
since it is observed primarily in the magnetospheric density cavity between the 
plasmapause and magnetopause (see an example in Figure 5). The trapped continuum 
spectrum is observed as a broad-band emission with very little frequency structure. The 
broad-band structure of the trapped continuum spectrum is believed to be produced from 
a series of narrow band emissions at slightly different frequencies from an extended 
source region at the plasmapause whose emission then mixes due to multiple reflections 
(with some Doppler broadening) in the magnetospheric density cavity. 

Continuum radiation at frequencies above the magnetopause plasma frequency has been 
referred to as the "escaping" component [Kurth et ah , 1981] since it propagates from the 
Earth’s plasmapause to well outside the magnetosphere (see an example in Figure 5). A 
common characteristic of all the escaping continuum radiation is that it has narrow 
frequency bands of emissions showing that the name continuum is not entirely 
descriptive of the radiation in this frequency range. Continuum radiation below 100 kHz 
is typically generated at the plasmapause in the pre-noon sector (06-12 LT) and has been 
called ct normal continuum” by a number of authors [e.g., Kasaba et al ., 1998; Hashimoto 
etal. y 1999]. 

Figure 6 has been adapted from Plate 3 and Figurel of Green and Boardsen [1999] and 
shows the observed emission cone of the normal continuum radiation by combining four 
years of Hawkeye plasma wave data overlaid with ray tracing calculations at that same 
frequency. These authors found that normal continuum radiation was generated in the 
pre-noon sector and is beamed into broad cones (half width of -35°) at low frequencies 
(-30 kHz). When modeled with ray tracing calculations, the resulting radiation pattern is 
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such that the trapped continuum radiation does not completely fill the magnetospheric 
cavity as previously believed. In addition, the trapped continuum component that consists 
of multiple reflections between the magnetopause and the plasmapause, as shown in 
Figure 6 (labeled as “trapped rays”), is confined to only -10° in latitude about the 
magnetic equator. 

Direction finding measurements from all four of the Cluster spacecraft from spin 
modulated Whisper data have confirmed many of the general properties of continuum 
radiation such as the plasmapause source regions [see Decreau et al ., 2004]. Near source 
region measurements have revealed new properties such as the importance of small-scale 
density irregularities in the local amplification of normal continuum and the possible role 
wave ducting may play in its generation. 

Normal continuum radiation at frequencies near the magnetopause plasma frequency has 
also been used to be an indicator of the global scale magnetosheath electron density. 
Nagano et al. [2003], using data from PWI on Geotail, observed a portion of the normal 
continuum spectrum that was trapped within the magnetsheath. The upper frequency was 
controlled by the f p in the solar wind after taking into account a time delay. The upper 
frequency was termed the transition frequency since it was the frequency at which the 
radiation changed from isotropic (due to multiple reflections in the magnetosheath 
randomizing its direction of arrival) to anisotropic (having a preferred direction). These 
authors found that a 0-3 hour time delay in the variation in the transition frequency of the 
normal continuum relative to the upstream f p may suggest that the magnetosheath plasma 
relaxes at the local Alfven speed rather than the convection speed of the solar wind. 

3.2 CONTINUUM ENHANCEMENT 

Large variations in continuum intensity over the 20 - 80 kHz frequency range lasting for 
up to 2 hours have been observed primarily in the midnight to dawn sector and were first 
reported by Gough [1982]. This distinct feature in the continuum emission spectrum has 
been called “continuum enhancement” by Kasaba et al. [1998]. The characteristic 
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emission spectrum of continuum enhancement is shown in the upper left spectrogram 
panel of Figure 5. Kasaba et al [1998], using simultaneous Geotail and Wind 
observations, showed that the normal continuum emissions in the pre-noon sector 
appeared to follow continuum enhancement events. This relationship has lead these 
authors to suggest that a series of storm time electron injections are generating both types 
of emissions. 

3.3 KILOMETRIC CONTINUUM 

Kilometric continuum (KC) is a major component of the escaping continuum radiation in 
the 100 - 800 kHz frequency range. Hashimoto et al [1999] have rediscovered this 
component of the emission and has sparked considerable interest in further understanding 
various aspects of this radiation that make it different from its lower frequency trapped 
and escaping (< 100 kHz) counterpart generated in the pre-noon sector and discussed in 
the previous section. The lower right hand panel of Figure 5 clearly shows the discrete 
emissions bands of KC extending from 17-24UT. The frequency range for kilometric 
continuum is approximately the frequency range of auroral kilometric radiation (AKR), 
but as shown in the lower right hand panel of Figure 5, there are significant differences 
that can be used to easily distinguish between these two emissions. KC has a narrow band 
structure over a number of discrete frequencies while AKR is observed to be a broad- 
band and sporadic emission and can be seen from 16-17:00 and from 21:30 to 24:00 UT 
in that spectrogram. 

Kilometric continuum has been observed at all local times, although it has been difficult 
to make a positive identification of the emission during the times when Geotail was in the 
late evening or early morning local time sector when AKR was active [Hashimoto et al., 
1999]. From Geotail and IMAGE observations Hashimoto et al. [1999] and Green et al. 
[2004a] have found that kilometric continuum is confined to a narrow latitude range of 
approximately []15 0 about the magnetic equator. Although these characteristics make it 
different from the lower frequency continuum discussed in Section 3.2, the similar 
spectral characteristics of the emission and its relationship to the plasmapause supports 
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the conclusion by Menietti et aL [2003] from Polar observations that the radiation is 
generated by the same mechanism. 

At lower frequencies, beaming of continuum radiation around the magnetic equator to 
latitudes as high as 50° has also been observed by Jones et aL [1987] (from 80 to 100 
kHz), Morgan and Gurnett [1991] (from 45 to 154 kHz), and by Green and Boardsen 
[1999] (from 24 to 56 kHz) as discussed above. The narrow beaming of kilometric 
continuum in magnetic latitude has made this emission difficult to observe routinely or 
for only short periods of time except for equatorial orbiting spacecraft with the proper 
instrumentation, such as Geotail. 

The source region for KC was originally identified by Carpenter et aL [2000] as coming 
from plasmaspheric cavities, but more recently Green et aL [2002] and Green et aL 
[2004a] clearly identified KC as being generated at the plasmapause, deep within notch 
structures that corrotate with the Earth. Figure 7 has been adapted from Figure 8 of 
Green et aL [2002] and Figure 1 of Green et aL [2004a] and illustrates the location of the 
KC source region and resulting emission cone pattern of the radiation to be consistent 
with the observations. Panel A of Figure 7 is a frequency-time spectrogram from the 
PWI instrument on GEOTAIL showing the banded structure of KC. Panel B shows the 
magnetic longitude versus the equatorial radial distance of the plasmapause (derived from 
the inserted EUV image of the plasmasphere) and the position of GEOTAIL during the 
KC observations of panel A. Panel C is a ray tracing analysis which shows that the 
structure of the plasmaspheric notch has a significant effect on the shape of the resulting 
emission cone through refraction. The process by which the notch structure is produced 
in the plasmasphere is not completely understood at this time. 

3.4 AURORAL KILOMETRIC RADIATION 

Over 35 papers have been published in the last 5 years on auroral kilometric radiation 
(AKR) demonstrating that after the classic paper by Gurnett [1974], AKR is, arguably, 
still the most intensely studied magnetospheric emission. It is beyond the scope of this 
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brief review to discuss all these results. Significant progress in understanding various 
aspects in the generation of AKR has been made in the areas of seasonal and solar cycle 
variations, polarization measurements, fine structure, and source region location. AKR 
research continues to be an important aspect of ionospheric-magnetospheric coupling 
during the substorm processes. 

The cyclotron maser instability [Wu and Lee , 1979] is believed to be the basic 
mechanism for the generation of AKR. For generation of AKR near the local electron 
gyroffequency, the cyclotron maser instability requires a source of free energy in the 
electron distribution and the following condition: 

f p /f g < 0.3 (1) 

The observation that the AKR emission frequency is very near the local f g in the auroral 
zone density cavity [e.g., Hilgers et al. , 1991] where the electron distribution function has 
a number of sources of free energy [Delory et al., 1998; Ergun et a/., 1998; Ergun et al, 
2000; Fung and Vinas , 1994] is the strongest observational evidence supporting the 
cyclotron maser instability. 

The AKR emission spectrum, in general terms, is observed over a broad frequency range 
from ~30 to -700 kHz with the peak intensity around 250 kHz. Kumamoto and Oya 
[1998] noted a seasonal difference in the AKR intensity from Akebono plasma wave data 
with increases in intensity in the winter polar region over the summer polar region. 
Kumamoto et al. [2003] expanded on these results and showed seasonal and solar cycle 
variations in the vertical distribution of the occurrence probability of AKR emissions 
from low altitude Akebono measurements. Using Polar and IMAGE wave data, Green et 
al [2004b] used the observed AKR spectrum to estimate the altitude range of the source 
region as suggested by Lee et al [1980] by mapping the observed frequency to the f g in 
the auroral density cavity. In both these studies, ionospheric density variations were 
believed to be mainly responsible for large-scale changes in the AKR source location. 
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Using two years of IMAGE RPI and one and a half years of Polar PWI data Green et al.- 
[2004b], found a dramatic frequency shift in the average AKR spectrum (intensity 
normalized to 8 Re) with dipole tilt. Panel A of Figure 8 clearly shows two main seasonal 
effects on the average AKR spectrum as observed by the IMAGE RPI instrument: 1) the 
AKR emission spectrum shifts down in frequency with increasing dipole tilt, and 2) the 
emission spectrum extends over a larger frequency range for negative dipole tilt angles 
than for positive dipole tilt angles. Panel B of Figure 8 is a qualitative illustration of the 
AKR spectrum, from panel A, mapped into contours of fp/f g for a 70° invariant auroral 
field line. Green et al. [2004b] found that the source region during extreme negative 
dipole tilt angles (winter) shows the auroral zone density cavity stretching from a region 
where f g = 500 kHz to nearly f g = 80 kHz. In contrast, for extreme positive dipole tilt 
angles the source region extent is much smaller, has moved further up the field line, and 
has a lower frequency boundary of 60 kHz and an upper frequency boundary near 250 
kHz. In addition. Green etal. [2004b] also found an overall solar cycle effect with the 
average AKR intensity lowered by as much as an order of magnitude during solar 
ma ximum . These results are consistent with Newell et al., [1996] who showed that there 
is a significant suppression of discrete aurora by sunlight and that ionospheric 
conductivity is a key factor in controlling the occurrence of discrete aurora. It was 
Gurnett [1974] who first clearly showed that AKR was closely correlated with discrete 
auroral arcs. 

The Polrad instrument on Interball 2 has made detailed polarization measurements of 
AKR distinguishing between L-0 and R-X polarization components. Figure 9, adapted 
from Hanasz et al. [2003], show a frequency time spectrogram of AKR from Polrad (top 
panel) and the polarization (middle panel) derived from determining the Stokes 
parameters [technique described by Hanasz et al 2000] and by factoring in the trajectory 
of the spacecraft relative to the source region (bottom panel). Hanasz , et al. [2003] has 
found that L-O polarized AKR, observed to be coming from dayside sources, is about 
three times more frequent than L-O polarized nightside sources which are observed about 
10% of the time. In addition, Hanasz et al. [2003] was able to distinguish two classes in 
the circular polarized AKR spectrum. The “regular” AKR is dominated by the R-X mode 
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and observed in the upper frequency portion of the spectrum and a weak L-O mode in the 
lower portion as shown in Figure 9. The “irregular” AKR is somewhat patchy with the 
two polarizations interwoven randomly over the whole frequency range. The reason for 
this variation in polarization is unknown. 

Another new spectral feature of AKR recently discovered by de Feraudy et al. [2001] 
and Hanasz et al. [2001] is an intense, nearly simultaneous broad frequency burst of 
AKR in as short as 6 s and as long as a few minutes that is generated during the fast 
expansion of the auroral bulge. These bursts of emission extend from 100 to 800 kHz 
and last for only a few tens of minutes. What remains to be explained is how the free 
energy for these waves be created so rapidly along such a large section of the source field 
line. 

The AKR emission spectrum also contains narrow band fine structure referred to as 
striations. Menietti et al. [2000] have analyzed a semi-random sample of the AKR 
spectrum using the high-resolution wideband instrument on the Polar satellite. New 
results from that survey show that the AKR fine structure, consisting of very narrow band 
emissions, is observed predominantly in the 40-215 kHz range and drifts in frequency 
with time. Negative drifting fine structures bands (decreasing frequency with time) are 
observed approximately 6% of the time, while the positive sloping (increasing in 
frequency with time) striations are observed at a much lower rate. The typical drift rate of 
the fine structure is observed to be between -8 kHz/s and -2 kHz/s and, for the positive 
sloping fine structure, between 5 or 6 kHz/s. In addition, there was a general decrease in 
drift rate with increasing frequency. Although the origin of the AKR fine structure is not 
known, it has been pointed out that the frequency dependence of the drift rate is 
consistent with AKR production stimulated by an upward propagating electromagnetic 
ion cyclotron wave [Menietti etal., 1996]. 

Applying a very long baseline interferometry (VLBI) technique to simultaneous 
observations from the four Cluster spacecraft, Mutel et al. [2003 and 2004] determined 
the source locations, at selected frequencies, of individual AKR bursts. The VLBI 
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technique involves precise timing of measured wave fronts from two or more locations 
separated by great distances. Six baselines can be formed from the four spacecraft 
yielding, when combined with the knowledge that the AKR emission is veiy close to the 
local gyrofrequency in the source region, source positions within an error of -200 km 
from individual burst events. From several thousand observed burst events, Mutel et al. 
[2003] found that the AKR sources were confined to the nightside auroral oval and that 
the distribution of AKR auroral source footpoints, during the same storm period (1-3 
hours long), had an overall spatial scale from 1000 to 2000 km. This is entirely consistent 
with a number of other earlier studies [e.g., Kurth et al., 1975; Gallagher and Gurnett, 

1 979; and Panchenko, 2003], As an example of the “pin-point” accuracy that this new 
technique has yielded. Figure 10 shows source locations from AKR burst at 125 kHz 
(red), 250 (green), and 500 kHz (blue) from two substorms. Panel A shows AKR burst 
locations during a southern hemisphere (positive dipole tilt) pass while panel B show the 
sources from a northern hemisphere (negative dipole tilt) pass. A clear shift in local time 
in AKR source region extent is seen in Figure 10 which is consistent with the results by 
Green et al. [2004b] that AKR sources are near dusk during times of positive dipole tilt 
and near midnigh t during times of negative dipole tilt. This effect is attributed to AKR 
source regions preferring the lower density cavity regions that develop along field lines 
whose foot is not illuminated by the sunlight. During times of negative dipole tilt a 
greater portion of the auroral oval is in darkness (AKR sources found near the dusk or 
late afternoon) with the opposite effect occurring during times of positive dipole tilt 
(AKR sources more towards midnight). It is important to point out that seasonal 
variations in the production of AKR were also found in an earlier study by Kasaba et al., 
[1997] from GEOTAIL/P WI observations. These authors showed that AKR is more 
active in the winter hemisphere for higher frequencies (500 kHz) representing lower 
latitude source regions. 

4. ACTIVE PROBES 

Plasma waves offer invaluable means for probing remote plasma regimes that are not 
readily or conveniently accessible to in situ measurements. Using the propagation 
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characteristics such as cutoffs and resonances of different wave modes [Stix, 1992], it is 
possible to ascertain plasma environment parameters, such as density and magnetic field 
strength profiles, with accuracies that rival or are even superior to those measured by in 
situ instrumentations. An advantage of active probing by plasma waves over conventional 
passive observation techniques is that a large volume of the target plasma regime can be 
probed simultaneously. In this section, we review recent advances in active remote 
sensing of magnetospheric plasmas. 

4.1 BRIEF HISTORICAL DEVELOPMENT 

The earliest applications of active probing of space plasmas were soundings of the 
ionosphere by radio waves [Appleton and Barnett , 1925; Breit and Tuve 1 926] . In fact, 
much of our current knowledge about the structures of the ionosphere has been derived 
from active radio measurements [e.g., Davies , 1990; Hunsucker , 1992]. Because of 
ionospheric shielding, ground-based radio observations of the magnetosphere are not 
possible except by whistler waves [ Storey , 1953; Helliwell , 1988.]! Nevertheless, active 
magnetospheric sounding has been proposed for many years [Franklin and Maclean , 

1 969; Ondoh et al , 1978; Calvert , 1981]. Because of the long ranges and uncertainties 
involving magnetopspheric plasma targets, the practicality of magnetospheric radio 
sounding had been subject of investigations [Green et al., 1993; Calvert et al ., 1995; 
Calvert et al, 1997; Greenwald, 1997a,b]. The first successful long-range 
magnetospheric sounder experiment is the Radio Plasma Imager (RPI) [Re ini sc h et al , 
2000] aboard the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) 
satellite [Burch, 2000], launched on March 25, 2000. 

4.2 ACTIVE RADIO TECHNIQUES 

Active probing by radio waves involves transmitting radio waves into a plasma target, 
and receiving the signals (echoes) returned after the incident waves interacted with the 
target. Classification of different probing techniques depends on the mechanisms by 
which signals are returned to the receivers, and thus on the frequencies of the probing 
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waves. The instrument that performs active radio probing of plasmas is called a sounder 
or radar. 

4.2.1 INCOHERENT SCATTERING 

As mentioned above, radio probing has been used most extensively in ionospheric 
studies. Notable instrumentations used in these studies include bottom-side (ionosondes) 
and topside (space-based) radio sounders, and incoherent and coherent scatter radars 
[Kelley, 1989]. An incoherent scatter radar (ISR) operating at high frequencies (> few 
hundred MHz) transmits waves with wavelengths that are typically small compared to the 
Debye length, h De =vJf v (where v e is the thermal speed and / p is the plasma frequency) of 
the intervening medium and the target plasma. The ISR-transmitted waves can therefore 
interact with individual charged particles and thus be Thomson scattered, i.e., absorbed 
and re-radiated. A small amount of the scattered wave energy can always return to the 
receiver as echoes. Signals from different ranges can be selected by different delay times. 
The ISR turns out to be a very useful tool for studying the ionosphere above the E and F 
layers where the densities are higher than the densities in the layers above. 

It should be noted, however, that when ISR was first constructed to observe the 
ionosphere, the detected signal strengths were found to be much stronger and 
concentrated in much narrower bandwidths than expected from incoherent backscattering 
[Bowles, 1958]. It was then found that the observed scattered signals at frequencies <few 
100 MHz were due to backscattering from ion acoustic waves, a process that is more akin 
to coherent scattering to be discussed next [e.g., Dougherty and Farley, I960]. 

4.2.2 COHERENT SCATTERING 

Ionospheric and magnetospheric plasmas are filled with irregularities of various scale 
sizes [e.g., Fejer and Kelley , 1980; Hanuise , 1983; Fung et al. , 2000; Fung et al . , 2003]. 
When the wavelengths of the probing radio waves are larger than the local Debye length 
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and become comparable to certain plasma irregularity scales, the waves can be 
effectively scattered by plasma collective effects according to 

kj = k s + ki,. (2) 

where ki, k s , and kir are the incident, scattered, and plasma irregularity wave vector, 
respectively. A special case of coherent scattering of this type is backscattering, in which 
k s « -ki, so that equation (2) gives k* ~ 2ki. This is the well-known aspect-sensitive 
backscattering condition that the scattering irregularity scale is half the probing 
wavelength. Coherent scatter radars (CSR) and aspect-sensitive scattering are particularly 
useful for measuring the occurrences of plasma instabilities and determining plasma 
irregularities, such as spread F [Fejer and Kelley , 1980; James , 1989]. 

4.2.3 RELAXATION SOUNDING 

Like any object, a plasma has a number of “ring tones” or resonances that are 
characteristic of the plasma density and embedded magnetic field strength. Plasma 
resonances have been explained in terms of sounder-stimulated electrostatic waves that 
propagate to short distances from the sounder and are reflected by local gradients back to 
the sounder receiver [ Muldrew , 1972]. To learn about the local properties of a plasma 
medium, a relaxation sounder can be used to send out small amount of wave energy to 
stimulate resonances in the plasma nearby to the sounder [e.g., Etcheto and Bloch , 1978; 
Etcheto et al, 1983; Trotignon et al. y 1986; Oya and Ono y 1987; Oya et al. y 1990]. This 
technique has been used by the Alouette-ISIS topside sounder satellites [ Jackson , 1986] 
to produce many valuable studies of ionospheric resonances [e.g., Benson , 1971; 1977; 
1996; Warnock , 1969; Warren and Hagg, 196%; Muldrew, 1972]. Typical plasma 
resonances observed in the ionosphere and magnetosphere include the electron 
gyroharmonic resonances n/ g , electron plasma resonance / p , and the upper hybrid 
resonance 


/ g (kHz) = 0.028 B(nT) 


( 3 ) 



6/22/04 


18 


/ p (kHz) ~ 8.97 [N e (cm' 3 )] 1/2 

(4) 

/uhr = (/g 2 +/ P 2 ) 1/2 

(5) 


As shown in Benson et al. [2003], gyroharmonic and other plasma resonances [(4) and 
(5)] are useful for deriving local magnetic field strengths and densities (to within a few 
percent), particularly when there is no other in situ observations of such parameters. 

In addition, there exist “(n+l/2]/g” resonances that occur between successive 
gyroharmonics and have been referred to as the Q„ and D n resonances, depending on 
whether they are observed above or below the upper hybrid resonance /„h [equation (5)] 
[e.g., Benson et al., 2003], Assuming a thermal electron plasma distribution, Warren and 
Hagg [1968] obtained an expression for the Q„ frequencies, 

/<> =/ g [n + (0.464/n 2 ) (/p// g ) 2 ] (6) 

For n > 1, the Q n resonances can be weak and not easily distinguishable from 
gyroharmonic resonances when (fp/f g ) 2 < 1 [Perraut et al., 1990]. Frequency splitting of 
Q n resonances has also been observed [e.g., Trotignon et al ., 1986; Canu et al., 2001; 
Trotignon et al ., 2001]. 

The identification and interpretation of magnetospheric D n resonances has been 
somewhat controversial [LeSager etal , 1998; Benson et al., 2001; Canu , 2001a,b]. D n 
resonances were first identified in topside ionospheric sounding data [Nelms and 
Lockwood , 1967; Oya, 1970], but there had not been any first-principle explanation for 
their mode of propagation. Osherovich [1987; 1989] and Osherovich and Benson [1991] 
have found, however, the eigen- frequencies derived from an electromagnetic force-free 
analysis of a cylindrical plasma volume to match the observed D n resonances quite 
successfully, yielding 
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/on = 0.95(/g/p) I/2 n 1/z 

(7) 

/Dn + =(&n 2 +/g 2 ) 1/2 

(8) 

s 

1 

CM 

£ 

It 

4 

(9) 


While equations (7)-(9) seem to allow an independent means to determine the local / g and 
f V9 particularly when the plasma and gyroharmonic resonances are not otherwise 
observed, their origin remain to be explained theoretically based on realistic plasma 
distributions. 


4.2.4 LONG-RANGE SOUNDING 

As implied by its name, long-range (LR) sounding uses radio waves to probe plasma 
regimes at a distance. This differs from the other techniques mentioned so far in that 
target ranges in LR sounding can be more than an Earth radius (1 Re ~ 6378 km). The 
long sounding ranges would normally require higher transmission powers to overcome 
the long distances. Special digital techniques, such as pulse compression and spectral 
integration, however, are used by RPI on IMAGE to lower the transmission power 
requirements so that LR sounding of the magnetosphere can become possible [Benson et 
al . , 1998; Reinisch etaL, 2000]. 

After a sounder pulse is transmitted, different types of echoes are received, depending on 
the structure of the medium and irregularities present in the target region [Fung et al . , 
2000], Figure 1 1 illustrates the different types of possible echoes from a magnetospheric 
sounder. As discussed above, short-range echoes from local resonances (e.g.,/ p and n/ g ) 
are due to electrostatic waves excited by the sounder pulse [Muldrew, 1972; Benson et 
a/., 2003 and references therein]. In LR magnetospheric sounding, due to the generally 
low transmission power involved, only specularly reflected or coherently backscattered 
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signals are detectable over long distances. Incoherent scatter signals are too weak to be 
observed. 

Radio waves transmitted at frequencies above the local cutoffs (f K or / P , see Figure 1) can 
propagate freely in all directions. When the wave propagates into a refractive index 
gradient (due to density or magnetic field increases), its wavelength will increase until it 
reaches a cutoff where the wave frequency matches the local cutoff frequency, and the 
refractive index vanishes (k = 0). The wave will then be reflected [Stfx, 1992]. According 
to the law of reflection and the wave-path reciprocity theorem [Born and Wolf j 1999], 
only waves propagating parallel to the refractive index gradient (k || Vn) will return to the 
sounder (receiver) location. 

If the target plasma has a smooth gradient with a scale length L^ B » X and is free of 
small-scale (A*. < A,) irregularities, then the incoming sounder waves will be specularly 
reflected at a cutoff (fR for R-X mode and f p for L-O mode as shown in Figure 1). The 
wave cutoffs are related to the local electron plasma density and magnetic field 
conditions, echoes of the sounder waves then yield important information on the remote 
target region. L-O mode echo frequencies yield remote density information directly. 

Since the R-X mode cutoff depends on both the electron density and magnetic field 
strength, estimates of remote densities can be obtained from R-X mode echoes if the 
magnetic field strengths at the remote sites can be determined from magnetic field 
models. The total time delay t of an echo then yields the so-called virtual range [R v = 
ct/2] of the reflection point, which assumes that the propagation speed throughout is the 
vacuum speed of light c. Techniques are available to invert the frequency- versus- virtual 
range observations into density profiles of the target plasma [e.g., Jackson , 1969; Huang 
andReinisch , 1982; Reinisch etal , 2001a]. 

On the other hand, if there are small-scale irregularities in the target region, then different 
frequency waves can be coherently backscattered wherever the ki r ~ 2k* condition is 
satisfied. The fact that the backscattering condition can be satisfied at different points in 
the target region means that there would be range spreading in the received echoes. In 
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addition, if the irregularities have both spatial and temporal variations (/ ir , k/), then 
analogous to equation (2) the scattered waves will also exhibit frequency spreading 
according to 

fi-fs+A ( 10 ) 

The amount of spreading will depend on the spectrum of the irregularities. Equations (2) 
and (10) represent the conservation of momentum and energy in the wave scattering 
process. 

It may be of interest to estimate the amount of frequency spreading due to electron 
plasma waves at the plasmapause where coherent aspect-sensitive backscattering are also 
observed [Carpenter et al y 2002; Fung et al, 2003]. Assuming a 200 kHz transmitted 
wave, the backscatter condition above then implies the presence of irregularities with a 
cross-field scale - 0.8 km. From the dispersion relations of electron plasma waves, 
either Langmuir waves or upper-hybrid waves, we have/ir 2 - f p 2 [ 1 +(3/2)(W^ r ) 2 ] or/u- 2 = 
f p 2 + f g \ where Xn = 7.43Te 1/2 N c ' 1/2 is the electron Debye length (in meter) with T e (eV) 
and N e (cm* 3 ) being the electron temperature and electron density, respectively. Thus for a 
T e ~1 eV and a N e -100 cm" 3 appropriate for the plasmaspause region, we have (A, D /ta) 2 
« 1 and > f p ~ 90 kHz, which implies that the scattered waves from the plasmapause 
can be changed by -50% of the transmitted frequencies. Given the prevalence of upper- 
hybrid waves in the plasmapause region [Carpenter et al.> 2002], we would expect to 
observe range spreading as well as frequency spreading in plasmapause echoes. 

The Radio Plasma Imager (RPI) on IMAGE performed the first successful LR 
magnetospheric sounding [Reinisch et al. y 2001b, Fung et al . , 2002]. The plasmagram in 
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Figure 12 shows a collection of echoes resulting from direct sounding by RPI near the 
equatorial plasmapause. The vertical bright lines are gyroharmonic resonances from 
which the local electron gyro frequency can be determined, / g = 16.34 kHz. Direct 
sounding of both the plasmapause and the outer plasmasphere resulted in spread echoes, 
indicating the presence of small-scale irregularities (A^ ~ 0.25 - 2 km) in these regions 
[Carpenter et al ., 2002; Fung et aL , 2003]. The discrete echo trace is due to signals 
guided along the geomagnetic field line and reflected from the high-latitude topside 
ionosphere. The similar amplitudes of the guided (discrete) and direct (spread) signals 
indicate that the guided signals have suffered very little path loss even though they have 
traveled much greater distances. 

Figure 13 depicts the signal guiding mechanism in which total internal reflections trap 
and help guide radio signals along the geomagnetic field in a region of enhanced index of 
refraction, much like a waveguide or optical fiber. Field-aligned irregularities (FAI) are 
especially efficient in producing guided echoes and have been observed in the 
ionosphere. [James, 1989; Hanuise, 1983]. Observations of guided echoes throughout the 
magnetopshere imply that FAI are widely present [Carpenter et al . , 2002; Fung et al., 
2003]. 

4.3 DETERMINATION OF MAGNETOSPHERIC PLASMA PARAMETERS 

Magnetospheric density and magnetic field strength can be determined by analyzing 
active sounding observations. We illustrate this by giving a few examples below. 

4.3.1 LOCAL DENSITY AND MAGNETIC FIELD STRENGTH 

Figure 14 shows observations of direct and guided echoes observed near the plasmapause 
(L = 3.8). The inset shows the IMAGE orbital position at the time of observation. The 
direct echoes were reflected from inside the plasmasphere along the direction 
perpendicular to the geomagnetic field lines. The guided echoes, on the other hand, 
propagated along the L=3.8 field line toward the high-latitude region and are reflected at 
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the /R-cutoff altitudes. It is important to note the nearly vertical <c up-side-down hook” 
feature of the guided echo trace. This feature indicates that the local /R-cutoff of the 
observation/transmission point was near 248 kHz, so waves at just above that frequency 
transmitted along the magnetic field propagate more slowly than waves at successively 
higher frequencies, leading to much longer virtual ranges (delay times) although the 
actual path is shorter. 

Based on the nominal local / g = 22.2 kHz determined from theTsygenanko-96 model 
magnetic field, the lines at 45, 66.6 and 1 10.9 kHz are consistent with being the 2 nd , 3 rd , 
and 5 th gyroharmonic resonances, which then yield an observed / g = 22.25 kHz. The 4 th 
harmonic was not detected because the programmed receiver frequencies (300 Hz 
bandwidth), each separated by a 4%-logarithmic frequency step, may not necessarily 
coincide with the harmonic frequencies. Using the /r and f g values, we can determine the 
electron plasma frequency to obtain /p = 236.6 kHz, consistent with the rather weak 
resonance observed at 233.6 kHz. The local plasma frequency thus determined 
corresponds to an electron density N e = 700 cm -3 . 

4.3.2 POLAR CAP DENSITY 

During IMAGE apogee passages over the northern polar region, polar cap echo traces 
have been observed to last up to three hours [Fung et al . , 2002; 2003; Nsumei et 
al., 2003]. In this high latitude region, the geomagnetic field lines are nearly radial, so that 
direct echoes from the polar cap and field-aligned guided echoes are not easily 
distinguishable. Inversion of polar cap traces (either from direct or guided echoes) will 
yield the nearly radial density profiles of the polar cap. Figure 15A taken from Nsumei et 
al. [2003] shows the radial electron density profiles of three successive soundings of the 
polar cap at 8 minutes apart. The three different starting high radial distances indicate the 
changing positions of the IMAGE spacecraft, reaching lower densities at the higher 
altitudes. The high frequency pulses at the same frequency apparently reached the same 
limiting altitude in the polar cap. 
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Combining successive observations taken over an interval, such as shown in Figure 15 A, 
it is possible to construct two-dimensional images of polar cap density in the IMAGE 
orbit plane, as shown in Figure 15B for two different geomagnetic (K p ) conditions. 

Figure 15B shows clearly the higher electron density content in the polar cap during 
geomagnetic active periods [Nsumei etal, 2003]. 

4.3.3 FIELD-ALIGNED DENSITY PROFILES AND STORM-TIME DYNAMICS 

It has long been recognized that the electron density distribution along the geomagnetic 
field line is an important parameter for understanding magnetospheric dynamics. Its true 
measurement has not been possible until magnetospheric radio sounding became 
successful [Reinisch et al, 2001a]. As described in Reinisch et al, [2001a] and also in 
Fung et al, [2003], guided or ducted echoes observed in the plasmasphere and 
plasmapause regions provide the only means by which field-aligned electron density 
distributions, like the polar cap case in Figure 15A, can be measured instantaneously. 
When guided echoes are observed to reflect in conjugate hemispheres, they sometimes 
form the epsilon signature [Fung et al, 2003], yielding the field-aligned density 
distribution between the northern and southern hemispheres. Figure 16A shows two field- 
aligned electron density profiles before and during a geomagnetic storm. The “quiet day 
model” was developed as a fit to derived field-aligned densities from a series of RPI 
plasmagrams, with ducted echoes, during one pass through the plasmasphere on the 
previous quiet day. The “measured” values are derived from one plasmagram of field- 
aligned echoes during the storm. Figure 6 A clearly shows the evacuation of plasma 
between ±40° magnetic latitudes at L=2.7 during the storm. 

During the course of a geomagnetic storm lasting 2-3 days, the 14.5-hr IMAGE orbit 
allows RPI to perform many soundings through the plasmasphere at different L values 
and to obtain the corresponding electron density distributions along the different field 
lines, revealing the associated plasmaspheric dynamics. Figure 16B shows the dramatic 
density decreases that occurred throughout the outer plasmasphere during the 
geomagnetic storm of March 3 1-April 2, 2001 [Reinisch et al , 2003]. The upper panel in 
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the figure shows the development of the storm as provided by the D st index throughout 
the pre-storm to post-storm interval. The lower panel shows the filling factor of different 
flux tubes (L) (normalized to their pre-storm densities) as measured by their 
corresponding guided echo traces, as a function of storm phases. The figure clearly shows 
the depletion or shrinkage of the entire outer plasmasphere during the stonn, with the 
formation of a new storm-time plasmapause near L pp -stonn = 2.5. Below this L value, the 
plasmasphere is largely unaffected, indicating the limit to the magnetospheric electric 
field penetration. The full plasmasphere is recovered at about 2.5 days after the main 
phase of the storm. 


5. CONCLUSIONS 

Within the last 5 years significant progress has been made in both understanding the 
location and importance of natural emissions in the inner magnetosphere and in using 
active EM waves as a method for remote sensing inner magnetospheric dynamics and 
structure. In this brief review several key advances have been discussed which include: 

1. The importance of lightning generated whistler mode waves in the plasmasphere has 
emerged as an important factor in maintaining the slot region in the electron radiation 
belts and becoming a component of plasmaspheric hiss. 

2. A more complex picture in the source location of non-thermal continuum radiation has 
emerged than previously thought. In addition to the normal trapped and escaping 
continuum radiation generated in the pre-noon section, continuum radiation at kilometric 
wavelengths is generated in plasmaspheric notches that corotate with the Earth, and 
continuum enhancement is observed in the pre-dawn sector. 

3. The AKR source region moves up auroral field lines and shifts toward midnight with 
increasing dipole tilt angle (summer) from a more low altitude and dusk location during 
the winter season. 

4. From direct and ducted magnetospheric sounder waves, the density structure and 
evolution of the plasmasphere and polar cap have been determined providing the 
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observations needed to understand the processes of plasmaspheric filling and the 
variability of polar cap outflows. 
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Figure Captions 


Figure 1. Trapped and escaping waves modes in a plasma in which the f p > f g illustrating 
the upper and lower frequency cutoffs and resonances. 

Figure 2. A summary of the four whistler mode waves in the inner magnetosphere as a 
function of their frequency range that are reviewed in this paper. 

Figure 3. Average plasmaspheric magnetic field spectral density showing frequency and 
radial distance distribution of whistler mode waves in the plasmasphere [after Andr€ et 
al., 2002]. 

Figure 4 Schematic of the longitudinal distribution of the whistler mode waves in the 
plasmasphere. 

Figure 5, The observed source locations of the escaping [after Dtfcrfau, et al. , 2004], 
trapped [after Gurnett and Frank , 1976], kilometric [after Hashimoto et al., 1999] and 
continuum enhancement [after Kasaba, et al., 1998] emissions. 

Figure 6, The observed emission cone of the trapped continuum radiation from four 
years of Hawkeye plasma wave observations at 3 1.1 kHz overlaid with ray tracing 
calculations at that same frequency [after Green et al., 1999]. The ray tracing calculations 
show the narrow range in wave normal angles about the magnetic equator in which the 
emission is reflected between the magnetopause and the plasmapause. 

Figure 7. The KC wave observations from Geotail/PWI (panel A) map to a 
plasmaspheric notch structure as observed by IMAGE/EUV (panel B) where the resulting 
emission cone pattern (panel C) is modeled with ray tracing calculations [after Green et 
al., 2002 and 2004a]. 
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Figure 8. The average spectrum of AKR from two years of data from IMAGE/RPI with 
respect to dipole tilt (panel A) clearly shows the dramatic frequency shift from high to 
low frequencies with increasing dipole tilt angle. Panel B illustrates the shift in the 
auroral zone density cavities corresponding to the observations in panel A [after Green et 
al, 2004b]. 

Figure 9. Detailed polarization measurement of AKR showing the L-0 and R-X 
polarization components [after Hanasz, et al, 2000]. 

Figure 10: AKR source region measurements from the Cluster/WBD instrument using 
all four spacecraft as a long baseline interferometer array. Panel A shows AKR burst 
locations during a southern hemisphere (winter) pass while panel B shows the shift in the 
AKR source measurements to earlier local times during a northern hemisphere (winter) 
pass [after Mutel, etal., 2004]. 

Figure 11: Different types of possible echoes resulting from direct long-range sounding. 
Guided echoes are produced by signals trapped by total internal reflection along field- 
aligned plasma irregularities. Direct (specular-reflected) echoes are produced by signals 
reflected by a wave cutoff, e.g at f=f v . 

Figure 12: Echoes resulting from direct radio sounding from outside the plasmapause (L 
= 5.0, 26° MLAT, 16.4 MLT), showing plasma resonances, scattered (spread), and 
guided (discrete) echoes [after Fung et al., 2003]. 

Figure 13: A sketch of field-aligned guided echoes trapped by total internal reflections. 
Field-aligned plasma irregularities create regions of enhanced refractive index n along 
geomagnetic field lines, causing radio signals to be guided by total internal reflections. 

Figure 14: Observations of local resonances and local cutoffs of guided and direct echo 
traces near the plasmapause allow the determination of the local plasma parameters / p and 

/ 8 - 
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Figure 15: (A) Polar cap radial density profiles obtained from successive RPI soundings; 
(B) 2-dimensional images of the polar cap structure during geomagnetic active and quiet 
periods [after Nsumei etal., 2003]. 

Figure 16: (A) Field-aligned electron density distributions at L= 2.8 before (the “quiet 
day model” is obtain from a fit of previous observations) and during a geomagnetic storm 
on April 1 , 200 1 ; and (B) the depletion of the entire outer plasmasphere (L > 2.5) and the 
establishment of a new plasmapauses during and after the geomagnetic storm [after 
Reinisch etal., 2003], 
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